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ABSTRACT 



We calculate the scattering of X-rays by interstellar dust, for a dust model that 

reproduces the observed wavelength-dependent extinction and polarization of starlight. 

<*' [ On interstellar sightlines that produce appreciable starlight polarization, we predict 

I that the dust-scattered X-ray halo around point sources will have measurable azimuthal 

' asymmetry due to scattering by partially- aligned nonspherical grains. We calculate the 
OO ' 

■ expected halo asymmetry. X-ray halo asymmetry provides a new test of interstellar 

^ ■ dust models. 

o 

Q^i Subject headings: ISM: dust, extinction - scattering - X-rays: ISM - X-rays: general 

o 
U 

^ ! 1. Introduction 



Interstellar dust grains scatter X-rays through small angles, as was first pointed out by Over- 
beck (1965), Slysh (1969), and Hayakawa (1970). Because of this scattering, an image of an X-ray 
! point source includes a "halo" of X-rays that have been scattered by dust grains near the line of 

sight. First observed by Catura (1983) using the Einstein observatory, scattered X-ray halos have 
since been measured by a number of telescopes, including Einstein (e.g., Mauche & Gorenstein 
1986), ROSAT (e.g., Predehl & Schmitt 1995), Chandra (e.g.. Smith, Edgar, & Shafer 2002), and 
XMM-Newton (e.g., Costantini, Freyberg, & Predehl 2005). 

The observed polarization of starlight (Hall 1949; Hiltner 1949) requires that interstellar grains 
be both appreciably nonspherical and partially aligned, with starlight propagating through the 
dusty interstellar medium becoming linearly polarized as the result of preferential attenuation of 
one of the linear polarization modes ("linear dichroism"). In addition to producing polarization 
of starlight, aligned dust grains produce polarized thermal emission at wavelengths from the far 
infrared to the microwave. Observers seeking to measure the polarization of the cosmic microwave 
background must subtract this polarized "galactic foreground" from observations of the microwave 
sky. Realistic models of nonspherical interstellar grains are therefore of interest for many reasons. 
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There continue to be many uncertainties concerning the composition and geometry of inter- 
stellar grains (for a recent review, see Draine 2003a), and X-ray absorption and scattering can be 
used to test grain models. We consider a specific grain model consisting of spherical carbonaceous 
grains and oblate spheroidal silicate grains, with size distributions and size-dependent degree of 
alignment adjusted to reproduce the observed wavelength dependence of both interstellar extinc- 
tion and interstellar polarization. 

The scattering "halo" produced by a nonspherical grain will not be azimuthally symmetric, and 
therefore the population of aligned interstellar grains can be expected to produce asymmetric X-ray 
scattering halos. Here we develop a method for calculating the differential scattering cross section 
for X-rays incident on grains with arbitrary geometry, and apply it to calculate X-ray scattering 
halos for a realistic model of interstellar dust. 

We find that the X-ray scattering halo produced by this model of aligned intersellar grains has 
appreciable and observable asymmetry. Wc propose statistics R^p to measure the asymmetry of 
observed halos. For two models of partially-aligned interstellar grains, we calculate the expected 
values of ¥6^"^ and , and discuss their observability. The predicted values of R^ should be 
measurable on sightlines where there is appreciable polarization of starlight. The ratio R^ /R^"^ 
is sensitive to the shape of the scattering grains. 

The paper is organized as follows: In §2 we discuss anomalous diffraction theory, as applied to 
interstellar grains, and our implementation of it. In §3 we discuss averaging over the grain rotation 
expected for partially-aligned suprathermally-rotating grains. In §4 we introduce statistics R^ to 
quantify azimuthal asymmetries in scattering halos. Wc obtain a realistic grain model in §5, and 
in §6 we use this model to predict R^ and R^ for sightlines where the magnetic field direction 
is uniform and perpendicular to the line-of-sight. We discuss using halo asymmetry measurements 
as a test of grain models in §7. Our results are summarized in §8. 

A reader concerned only with the predicted observability of this phenomenon may choose to 
proceed directly to §7. 



2. X-Ray Scattering by Dust: Anomalous Diffraction Theory 

2.1. Defining the Scattering Problem 

We consider an incident monochromatic electromagnetic plane wave with time dependence 
g-iu)t ^j^Q^ explicitly mentioned in the subsequent formulation) and spatial dependence exp(ik • r) 
which is fully coherent over the grain volume and totally non-coherent over the length-scales of 
the distances between them. The grain is assumed to consist of material characterized by complex 
refractive index m. We neglect effects arising from the crystalline or amorphous atomic structure 
of the grain (Bragg diffraction, for example, affects the scattering only at large angles, where the 
scattering halo is impossible to observe). 
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An electromagnetic wave, in a homogeneous (and nonmagnetic) medium, would propagate as 



As usual, the properties of the vector field can be understood by studying the behavior of each of 
its components, hence reducing our analysis to a scalar theory of scattering, governed by a scalar 
wave equation V^?7(r, Lo) + k'^m'^U (r, lo) = 0. The difficulties arise from the coupling of the different 
vector components by the boundary conditions at interfaces between different media. 

The scattering problem hence reduces to finding how the scatterer responds to an incoming 
wave [Zinc, producing an outgoing wave J/gcai where U = CAnc+t^sca should satisfy the wave equations 
above, as well as appropriate boundary conditons. For spheres, an exact series solution, first 
described by Mie (1908) and Debye (1909) and commonly referred to as "Mie theory", can be 
employed provided the sphere is not too large relative to the wavelength of the incident radiation. 
However, nonspherical targets require other methods. 



At X-ray energies materials have refractive indices very close to unity (|m — 1| <C 1) and the 
dust grains responsible for most of the scattering are usually much larger than the wavelength of 
the incoming radiation {ka » 1). In this regime, the scattering and absorption of X-rays can be 
calculated using an approximation first developed by van de Hulst (1957) and known as "anomalous 
diffraction theory", hereafter ADT. 

ADT is a combination of ray-tracing optics (applicable because ka S> 1) and Huygens' principle 
of propagation of a scalar field, applied to cases where the EM wave can enter and propagate through 
the grain with essentially no reflection or refraction (|m — 1| ^ 1). 

Consider an incident plane wave Uinc = Uqc^^^, propagating in the z direction. Under the 
conditions above, the plane wave, once reaching a plane V located just beyond the grain and 
normal to the direction of the propagation of the incident wave, will have changed by a fractional 
amount which we refer to as the shadow function, f{x,y): 



V^E(r,a;) + A;^m^E(r,u;) = 



(1) 



2.2. Anomalous Diffraction Theory 



f{x,y) = 1 -exp[i$(x,y)] 



(2) 



where the complex phase function $ is 




(3) 



where m{x, y, z 



) is the refractive index at the point {x,y,z). 



Once the shadow function is known, Huygens' principle allows the amplitude of the scattered 
part of the wave, in the radiation zone, to be calculated as a Fourier transform of the shadow 
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function over the plane V: 

UUrn) = C/o^^^^5(n) , (4) 
k"^ f 

S{n) = S{kx,ky) = — cxp[i{kxX + kyy)]f{x,y)dxdy , (5) 

2tt J 

kx = A;(n-x) , (6) 
ky = k{h-y) . (7) 

The scattering properties of the grain can be obtained from S{n), with the differential scat- 
tering cross section given by: 

The extinction cross section can be obtained from the optical theorem [see, e.g., Bohren & Huffman 
(1983)] 

(Text = ^MSm (9) 

= (l-e-*2cos$i) dxdy (10) 

where $i = Re($), $2 = Im($). Radiation traversing the grain at {x,y) is attenuated by a factor 
g-2*2^ so that a fraction (1 — e"^*^) of the incident power/area at {x,y) is absorbed by the grain. 
Thus the absorption cross section in the geometric optics approximation assuming |m — 1| <C 1 
(reflection and refraction at interfaces is small) is 

(Tabs = y (1 - e-'*^) dx dy , (11) 
and the total scattering cross section o-gca = o"ext — o"abs is 

(Tsca = J [l-2cos$ie-*2 +e-2*2] dxdy (12) 

= j dx dy . (13) 

The dimensionless efficiency factors for absorption, scattering, and extinction arc here defined to 
be Qabs = o-abs/Tra^g, Qsca = o-sca/vra^jj, Qcxt = Qabs + Qsca = o^ext /TTa^jj , where for a grain of solid 
volume V, the effective radius Oeg = (3y/47r)^/^ is the radius of an equal- volume sphere. 



2.3. ADT for a Sphere 

For a sphere and scattering angle 6, the simple form of the shadow function allows S{Q) to 
be written (van de Hulst 1957) 

S{e) = {kaf du(l-e-'''''"'^) Jo(/cae cos «) Sinn cos u , (14) 
Jo 
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p = 2ka{m-l) , (15) 

where Jq is the Bessel function of order 0. The extinction, absorption, and scattering efficiency 
factors are 

^ , 4 {cos 2/3 - e-P^ [cos(pi - 2/?) - \p\ smjp, - (3)]} 

yext = 2+ |-p (Ibj 



e 



-2p2 f>-2^)2 _ 2_ 



Qabs = 1 + + ^r^ — (17) 

P2 2p^ 

Qsca ~ Qext Qabs (1^) 

where pi = Re(p), p2 = Iin(p), and /3 = arctan(p2/pi). 



2.4. The DADT code 
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Fig. 1. — Differential scattering cross section for MgFeSi04 silicate spheres with radii o = 0.1 and 0.2/Ltm at 
E — 0.5, 1 and 2koV calculated with Mie theory (solid curve) and DADT (dots). The two methods coincide to within 
the accuracy of the plot. The refractive index sX E = 0.5, 1 and 2keV is taken to be m = 1— 2.079 x 10~^+3.201x 10~^i, 
1 - 7.152 X 10"* + 1.887 x lO'^i, and 1 - 1.920 x 10"" + 2.807 x lO'^i, respectively (Draine 2003b). 
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For general shapes, eq. (5) requires finding f{x,y) numerically, followed by a two dimensional 

integration for each scattering direction n. Because we will typically be interested in many scat- 
tering directions, it is advantageous to employ fast Fourier transform (FFT) methods to find S'(n). 

Our ADT-codc, hereafter referred to as Discrete-ADT code, DADT for short, samples the shadow 
function on a (x, y)-grid of 2^^ x 2^^ points, with 2"^ x 2'^ points within the projected area of 
the grain. ^ A 2-dimensional FFT then yields S{kx, ky) over a rectangular (kx, ky) grid of 2^^ x 2^^ 
points. The single-precision FFT code GPFA, developed by Temperton (1983, 1992) is employed 
for the 2-D FFT. The resulting S{kx, ky) is then transformed onto a 2^^ x 2^^ polar lattice (equal 
number of divisions in (j) and 0) using a two-dimensional cubic spline. 

We have tested DADT by comparing the scattering halo pattern calculated using DADT to results 
calculated with the Mie theory implementation MIEVO by Wiscombe (1980, 1996), after conversion 
to double precision arithmetic. Figure 1 shows the differential scattering cross section calculated 
for a = 0.1 and 0.2;[im silicate spheres at E = 0.5, 1, and 2kcV. The results calculated with DADT 
and those calculated with MIEVO are indistinguishable. Because MIEVO and DADT follow entirely 
different approaches to the calculation, this confirms the accuracy of both for the cases considered. 
Note that absorption is strong in some of the cases shown: Im(m)A;a = 1.62 for a = 0.2/xm and 
E = 0.5 keV. ADT requires |m — 1| <^ 1, but Im(m)A;a need not be small. The validity condition 
\m — 1 < 0.01 is fulfilled for sihcates for E > 250 eV. The ray optics validity condition ka > 10^ is 
satisfied for a > 0.08/xm(250eV/£;). 

2.5. Results 

Having verified that DADT yields accurate results for X-ray scattering by spheres, we now apply 
it to calculate X-ray scattering by oblate spheroids. Figure 2 shows the differential scattering cross 
section for axial ratio b/a = \/2 oblate silicate spheroids with symmetry axis perpendicular to 
the line-of-sight for two sizes {ues = 0.1 and 0.2/xm) and three energies {E = 0.5, l,2keV). The 
differential scattering cross section da/dQ,{@, 0) is shown as a function of 9 for = (the direction 
of the short axis of the grain) and cf) = 90°. It is apparent that the scattering is more extended 
in the (p = direction: it is easily shown that for this orientation of an oblate spheroid, the halo 
extent in the ^ = direction is larger than the extent in the cp = 90° direction by exactly the factor 
b/a. 

Figure 2 also shows that as a^s or E are increased, the angular extent of the halo shrinks, 
scaling as l/(aeff-E). 



^The shadow function f{x, y) is nonzero only over the projection of the target onto the {x, j/)-plane. However, any 
discrete Fourier transform requires that the range of integration be extended to distances r^ato evaluate ^(n) for 
angles within the first minimum of the scattering halo. These angles play a crucial role in our understanding of the 
halo properties, and we therefore extend our grid by a factor 2* in each direction, providing an angular resolution of 
~ 1/10 the angle of the first minimum of >S'(n). 
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Fig. 2. — Differential scattering cross section for h/a = \f2 oblate silicate spheroid with its symmetry- axis perpen- 
dicular to the line of sight and in the direction = 0, for aeff = 0.1/xm and 0.2/im. Results are shown for = 0° 
(solid line) and = 90° (dotted line). Scattering is more extended in the direction parallel to the "short axis" of the 
shadow function. 



Figure 3 shows the scattering for the same two grains as in Figure 2, but now with the grains 
oriented with symmetry axis parallel to the line-of-sight. For this orientation the scattering is 
azimuthally symmetric. The angular extent (e.g., locations of minima and maxima) is similar to 
the scattering in the ^ = 90° direction in Fig. 2, being proportional to 1/6 in both cases. 

Figure 4 shows contours of constant da/cKl for oblate spheroids, with symmetry axis (short 
axis) in the x direction. The contours are ellipses with axial ratios h/a. 



3. Scattering by Spinning, Precessing Grains 

In this section we give a few general results concerning symmetries in scattering by populations 
of spinning grains precessing around the local magnetic field. Rotation and precession are assumed 
to be rapid enough to ensure phase averaging over rotation angle and precession angle. 
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Fig. 3. — Same as Fig. 2, but for the grain symmetry axis parallel to the line of sight. Results are shown for = 0° 
(solid line) and 4> = 90° (dotted line); the two curves coincide because of symmetry. 

Consider radiation propagating in the z direction. Suppose that there is a magnetic field B 
in the x — z plane; let 9bz be the angle between B and z. The scattering direction is defined by 
scattering angles (6,0), where is the deflection angle, and cf) is an azimuthal angle with cf) = 
corresponding to scattering in the ic — z plane. 

Scattering by a sphere, or by a grain that is rotationally symmetric about the line of sight, will 
be azimuthally symmetric. This will also be the case for scattering by a population of arbitrarily- 
shaped grains if the distribution function for grain orientations is azimuthally symmetric about the 
line of sight. 

Nonspherical grains that are not randomly oriented will produce scattering that will not be 
azimuthally symmetric. Grains will in general be spinning rapidly; let J be the instantaneous 
direction of the grain angular momentum. The torques acting on a spinning grain have been 
discussed elsewhere (e.g., Draine Sz Weingartner 1997). The spinning grain will acquire a magnetic 
moment n \\ —J due to the Barnett effect. The /x x B torque drives precession of J around B. 
The precession period is short, of order weeks, so that we may assume a uniform distribution of J 
around the precession cone. The distribution of directions J will therefore be determined by 0bz 
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Fig. 4. — Contours of constant [dc7/df2(e, (/>)]/[dCT/df2(0, 0)] = 10""/^ for grains with b/a = ^2 (left panel) and 
h/a = 2 (right panel), with short axis in the x direction. Contours at local maxima are labelled; regions where the 
contours arc closcly-spaccd arc minima. The contours should ideally be perfect ellipses; deviations result from the 
discreteness of our grid and inaccuracies of interpolations by the contour-plotting software. Results were calculated 
for Oeff = 0.2/ttm silicate grains and E = IkeV, but the plots apply to other values of Ues, E, and composition 
provided the ADT validity criteria |m — 1| <C 1 and ka^ 1 are satisfied and, in addition, \m — l\ka <C 1. 



and the angle Gbj between B and J. 

We will assume the aligned grains in the present study to be spinning suprathermally, i.e., with 
rotational kinetic energy E'rot ^ ^Tg^, where Tgr is the grain temperature. For a suprathermally- 
rotating grain with fixed angular momentum J, dissipation resulting from viscoelasticity, the Bar- 
nett effect, or nuclear magnetism (Purcell 1979; Lazarian k, Draine 1999a,b) will cause the grain to 
minimize its rotational kinetic energy, and therefore to be oriented with ai || J or ai || —J, where 
ai is the principal axis of largest moment of inertia. The states with ai || J and ai || —J have 
the same energy, and we will assume them to be equally occupied. Therefore, for a suprathermally 
rotating grain, the distribution of ai in space will be fully determined by 0bz and the alignment 
angle 6b j- 

It is convenient to define a normalized scattering function 

^{QsAs) = -^(%^) (19) 

a,ca = /^^(^) ' (20) 

where (Gs,^s) are scattering angles, {dasca,/d^) is the differential scattering cross section per H 
nucleon summed over the different grain types and averaged over the ensemble of grain orientations. 
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and (Tsca is the total scattering cross section per H nucleon. With this definition, J ad^ = 1. For 

suprathermally-rotating grains, the normahzed scattering function a{Q,(j)) is determined by the 
grain properties (composition, size, and geometry), the angle 0bz and the distribution function for 
the alignment angle Qbj for each grain type. 

The function a will satisfy the symmetry^ 

a(e,0) = a(e,-</.) (21) 

and can therefore be written 

CO 

aiQ,^) = ao(e) + 2^af(e)cos(^(^) , (22) 

aeie) = — #cr(e,(/.)cos(£(/.) . (23) 
27r Jo 

If ©Bz = or 7r/2, the scattering will be symmetric upon reflection through the y axis: 

a(e,(/.) = a(e,7r-<^) , (24) 
ai = for odd ^. (25) 

For < 0Bz < 7i"/2, the symmetry (24) is not strictly required for scattering by dust. However, 
when anomalous diffraction theory applies, the symmetry condition (24) will also apply,^ in which 
case cr(0, 0) can be written 

oo 

a{e,(f>)=ao{@) + 2Y,MQ)cos{2e(t)) . (26) 

i=i 



4. Azimuthal Asymmetry in X-Ray Scattering by a Nonspherical Grain 

4.1. Asymmetry of Grain Scattering 

Consider some specific magnetic field direction Qbz, and an ensemble of grains with some 
specified distribution function for the alignment angle Gb j . The function 

/'02 /"StT 

g^''\ei,@2) = / d0sine / d4> a {&,(!)) (27) 
Jei Jo 



^Eq. (21) need not be satisfied for arbitrary grains, but will be satisfied if every grain in the population has a 
mirror-image counterpart. 

^This is a consequence of the scattering depending only on the shadow function (2). For any J, rotation of the 
grain around J and equal numbers of grains with ai || J and ai || —J together ensure that the ensemble of shadow 
functions is symmetric under reflection x — > —x, in which case a{Q,4>) = a{Q,TT — 4>). 
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Fig. 5. — M^^O, e) and R'f\o, &) vs. QaenE for b/a = V2 and b/a = 2 silicate spheroids. Results are shown for 
tteff = 0.1 and 0.2/im, and for E = 0.5, 1, 2keV. 



2tt I desmeao(e) 
101 



(28) 



gives the fraction of the scattering that is within the annulus [0i, 02]- 

For scattering angles in the annulus [6i, G2] the degree of azimuthal asymmetry of the differ- 
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B. T. Draine 2006.07.08. 

Fig. 6. — Geometry for scattering of X-rays by a dust grain. Angles are exaggerated; actual halo angles 9 are < 1°. 



ential scattering cross section can be characterized by functions 

pW,^ _ Jg^^e sin 9 cosm a{e,<f>) 

Jq^ ciB sm Jq c(0 (p) 

f®' dGsinG a^e) 
= — - for£>l . (30) 

/g'desine ao(e) 

The symmetry condition (24) yields R^^"^ = for odd i. 

We will be interested primarily in the quadrupolar asymmetry, characterized by the function 
R^^ ■ Figure 5 shows R^\o,Q) for perfectly aligned oblate spheroids. r'^\o,@) rises rapidly 
from zero to a value ~ 0.17 for b/a = \f2 and ^ 0.33 for hja = 2. We also see that the oc- 
tupolc/quadrupole ratio R^^^ / R^^ is sensitive to the grain shape: for Qa^sE > 3.5 arcmin/xmkeV, 
the ratio R^^ /R^^ increases from 0.17 to 0.33 as b/a is increased from ^/2 to 2. 



4.2. Asymmetry of X-Ray Scattering Halo 

The geometry of X-ray scattering is illustrated in Figure 6; note that X-ray scattering is 
significant only for small scattering angles 9, and therefore only for small halo angles 6. If the 
dust grain distribution in space is p{z) (assumed, for the moment, to be independent of x and y for 
small displacements from the linc-of-sight), then, for a steady isotropic source with luminosity per 
unit frequency L^,, the specific intensity of singly-scattered photons is given by (see, e.g., Draine &; 
Tan 2000) 
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p{z = CD) 



eg{(,9) = 6' + arctan 



(1-C)tan0 



C 



Tsca = -/VrCTs, 



(32) 

(33) 
(34) 



provided Tgca ^ 0.3 so that multiple scattering can be neglected. Let I{9, (j)) be the observed 
intensity of scattered X-rays. The function 



27r 



J^'UesmO J^"" d(l) I{9, 



(35) 



gives the fraction of the total scattered power that falls in the annulus [^1,^2]- The azimuthal 
asymmetry of the scattered halo in the annulus [^1, 62] can be measured by the function 



uiDra a ^ - fel dO siuO j^'' dct> cos(^</.) 7(0, <^) 



r-27r 



j;^d9 Sine j^''d4>m 



2tt 



(36) 



For optically-thin scattering by dust aligned by a uniform magnetic field perpendicular to the 
line-of-sight, we can calculate Rf^- 



J,'^detan9jidC aejOs) /5(C)/[C^ + (1 - Q^tan^^] 
detune dC aoiSs) piO/iC + (1 - C)2tan2 9] ' 



(37) 



where the are given by eq. (23). If we assume p{Q = ioi ( <^ 1 (i.e., negligible scattering from 
dust very near the source) we can assume tan^ <C 1, ^ ~ C^sj and approximate (37) by 



li dC m le^/c aoie) 9 de 



(38) 



5. Models for Aligned Interstellar Grains 

5.1. Observed Polarization 

Our objective is to calculate X-ray scattering from a realistic model of partially aligned non- 
spherical dust grains, with the size distribution and degree of alignment constrained to reproduce 
both extinction and polarization as a function of wavelength. 

The observed polarization of starlight (Hall 1949; Hiltner 1949) demonstrates that interstellar 
grains are both nonspherical and systematically aligned. We take the polarization as a function of 
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Fig. 7. — Extinction for randomly-oriented grains with the size distribution shown in Figure 9. Dotted curve shows 
the observed extinction law that was applied as a constraint on the model. The extinction for the models with 
b/a — \/2 (solid line) and 2 (broken line) nearly coincide and cannot be distinguished in this plot. 

wavelength to be empirically described by the "Serkowski law" (Serkowski 1973) for A < A^;, and a 
power-law for A^; < A < 5/im: 



The wavelength of peak polarization has a typical value Amax ~ 0.55;um, but varies from one 
sightline to another. The value of K is correlated with Amax) with K ~ 0.91(Ainax/0-55/um) + 0.01 
(Whittet et al. 1992). We take /3 = 1.7, within the range observed by Martin et al. (1992). For 
Amax = 0.55/im, we have K = 0.92, and A^. = 1.39^m. 
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Fig. 8. — Polarization cross section (Tpoi(A) for a line of sight perpendicular to the local magnetic field, for oblate 
silicate spheres with the size distribution of Fig. 9 and size-dependent alignment fraction of Fig. 10. Dotted curve is 
the observed polarization (eq. 39,40) that was applied as a constraint. The models with hja — \/2 (solid curve) and 
2 (broken curve) nearly coincide in this plot. 



Studies of many sightlines (Serkowski et al. 1975) find that 

< 0.09£;(S - y)/mag 



(42) 



or 



Pmax < 0.028y4(Amax)/niag . 



(43) 



Sightlines with Pmax/^(Amax) < .028/mag are assumed to pass through regions where the magnetic 
field direction is not transverse to the line-of-sight, or varies along the sightline, or where the degree 
of grain alignment is for some reason lower than average. 
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5.2. Models With Partially-Aligned Dust Grains 

We model the dust as a mixture of carbonaceous particles (including PAHs) and amorphous 
silicate particles. Models of the infrared emission from interstellar dust (Li & Draine 2001) require 
a population of PAH particles containing CpAH/Htotai ~ 30 — 60ppm. The present models assume 
PAHs to be present with a CpAn/Htotai ~ 55ppm, with optical properties as described by Li 
&: Draine (2001). In order to reproduce the observed wavelength-dependent extinction, larger 
grains composed of both amorphous silicates and carbonaceous materials are required (e.g., Mathis, 
Rumpl, & Nordsieck 1977; Draine & Lee 1984; Weingartner & Draine 2001; Zubko, Dwek, k. 
Arendt 2004). The carbonaceous particles are here taken to be spherical or randomly-oriented (and 
therefore not contributing to polarization), and the silicate particles are assumed to be spheroids, 
with diameter 2a along the symmetry axis, and 2h perpendicular to the symmetry axis. The size 
of a spheroid will be characterized by a^s = (a6^)^/^. 

For the spherical carbonaceous grains, the IR to UV extinction was calculated using Mie theory, 
using the dielectric tensor of graphite with the usual 1/3-2/3 approximation (Draine & Malhotra 
1993). For the silicate grains, taken to be oblate spheroids, we use the extended boundary condition 
method (EBCM) introduced by Waterman (1971) and developed by Mishchenko & Travis (1994) 
and Wielaard et al (1997) - see the general review by Mishchenko, Travis &; Mackowski (1996). 
Our computations make use of the code ampld.lp.f (Mishchenko 2000).^ EBCM codes encounter 
computational difficulties when the target becomes large compared to the wavelength. For the 
silicate oblate spheroids considered here, the ampld.lp.f code appeared to converge for b/X < 3.88, 
but sometimes failed for larger values of 6/A. Thus for b/X < 3.88 we used ampld.lp.f to calculate 
Qext for different orientations, but took 

Qext{b/X, b/a, m) Qext(3.88, b/a, m) for b/X > 3.88 . (44) 

For large values of b/X, and for refractive indices appropriate to the optical and ultraviolet, the 
extinction tends to be close to twice the projected geometric cross section (with zero contribution 
to polarization of starlight in this limit). For the size distributions characteristic of interstellar 

grains, grains with b/X > 3.88 make only a minor contribution to the total extinction, and a very 
small contribution to the polarization of starlight, so the approximation (44) does not introduce 
significant error. 

We consider models where all the silicate particles are oblate^ (6 > a) spheroids with a single 
axial ratio b/a, independent of size. The degree of alignment of grains of size a^s is given by an 
alignment fraction /(aegf ); < / < 1, where / = for random alignment, and / = 1 for grains where 



^ampld.lp.f is available at http://www.giss.nasa.gov/~crmim/t_matrix.html 

'There is some indication that oblate shapes provide a better match to polarization observations (Draine & Lee 

1984). Oblate spheroids also have the advantage of being invariant under rotation about their principal axis of largest 
moment of inertia, thus eliminating the need for averaging over grain rotation. 
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Fig. 9. — Size distributions for carbonaceous spheres and amorphous sihcate oblate spheroids. Size distributions are 
constrained to be smooth, and to reproduce the observed extinction and polarization. The size distributions shown 
employ about twice as much silicate material as would be permitted by current estimates of solar abundances. The 
size distributions are not shown for a^s < 0.01/im, because those grains contribute negligibly to X-ray scattering. 

the short axis (ai) is perfectly-ahgned with the magnetic field direction. We assume the magnetic 
field to be perpendicular to the line-of-sight, and for simplicity we approximate the distribution of 
partially-aligned grains using "picket fence alignment" : if the line-of-sight is in the z direction, and 
the magnetic field is in the x direction, we assume that a fraction (1 -|-2/)/3 of the oblate spheroids 
have ai || x, (1 — /)/3 have ai || y, and (1 — /)/3 have ai || z. This is not a physically realistic 
distribution - as discussed in §3, one should properly integrate over some distribution function 
(/)(0Bj) for the angle between the grain angular momentum J and the magnetic field direction B. 
However, given our lack of knowledge of the functional form of (pi^Bj), it is reasonable to assume 
picket-fence alignment, with a single number /(a) characterizing the alignment of grains of size a. 

We carry out a nonlinear least-squares fit for three continuous functions: the size distributions 
{dn/da)ca,T and {dn/da)^\\ for carbonaceous grains and silicate grains, and the alignment function 
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Fig. 10. — Fractional alignment for oblate silicate spheroids, as a function of effective radius a^g. Large grains are 
perfectly aligned, with their principal axis of largest moment of inertia parallel to the local magnetic field direction. 
Small grains are only minimally aligned. 



/(ocff). This fit is carried out using a number of constraints, embodied in a penalty function P = 
that we seek to minimize. One of the constraints is the requirement that iJmax/^(Amax) = 
.028/mag, with p{\)/pma.x required to reproduce eq. (39,40). 

In addition to a penalty for deviations from the observed extinction and polarization (at 100 
wavelengths, logarithmically spaced from A = 2.5/im to A = 0.111/xm), the penalty function P 
includes terms that are designed to favor solutions for which the size distributions dric/da and the 
alignment function /(a) are smooth functions of a. In addition, the penalty function includes a 
term designed to favor alignment functions /(a) with df /da > 0, as it is expected that small grains 
will be disaligned because of the effects of "thermal flipping" (Lazarian & Draine 1999a), while 
starlight torques (Draine & Weingartner 1997) will allow large grains to overcome the effects of 
thermal flipping and achieve suprathermal rotation, resulting in efficient alignment of large grains 
with the interstellar magnetic field in diffuse regions. 



-19- 



In addition, we include a penalty if the size distribution uses more than the amount of silicate 
and carbon that would be present in solid form in the ISM if the total interstellar abundances of C, 
Mg, Si, and Fe were equal to current estimates of solar abundances. Because we have no accepted 
theoretical expectation for either the size distributions dn/da or the alignment function /(a), the 
form of the penalty terms, and the weights they are given, are necessarily somewhat arbitrary. Our 
adopted penalty function P is described in the Appendix. 

Kim k. Martin (1995) modeled the polarization of starlight as a function of wavelength using 
silicate spheroids. For oblate spheroids with axial ratio h/a = \/2, they found that the polarization 
can be explained if the grains with a^s ^ 0.1/xm are nearly perfectly aligned, in which case a 
sightline perpendicular to B will have polarization Pmax/^y ~ 0.028/mag, the maximum value of 
the ratio observed in the interstellar medium. For larger values of 6/a, only partial alignment of 
the a > 0.1/um silicate grains is required to produce the observed polarization. 

Following Kim k, Martin wc take the silicate grains to be oblate spheroids with axial ratio 
b/a = -v/2, and we find the size distributions {dn/da)cax and {dn/da)si\ and alignment function 
/(«efr) giving the best agreement with the observational constraints. We repeat this exercise for 
silicate grains with b/a = 2. 

Figure 7 shows the extinction calculated for the two dust mixtures, together with the observed 
extinction. Figure 8 shows the polarization as a function of wavelength for the model, as well as the 
observed polarization. Both models successfully reproduce the observed extinction and polarization. 
However, in order to do so, both models require about twice as much silicate material as would 
be allowed by current estimates for the solar abundances of Mg, Si, and Fc. Sofia & Meyer (2001) 
recently discussed the applicability of solar abundances to the interstellar medium, and argue that 
abundances in young F and G stars provide a better standard, with Mg/H and Si/H about 12% and 
23% above the values currently favored for the Sun (although the Mg and Si abundances in F and 
G stars have large uncertainties - 40% and 33%, respectively), but this is still less than the amount 
of Mg and Si required to reproduce the observed extinction and polarization for the grain model 
considered here. This abundance shortfall is generally encountered by models that use "compact" 
grains to reproduce interstellar extinction (e.g., Weingartner & Draine 2001). It may indicate that 
the abundances of interstellar Mg and Si are higher than estimated from either the Sun or young 
F and G stars. Alternatively, it may indicate that interstellar grains have other geometries, e.g., 
"composite" grains with vacuum fractions of order 50% or more (e.g., Zubko et al. 2004), which 
might permit the observed extinction to be accounted for using less material in grains. 

In Figure 9 we show the best-fit size distributions of silicate and carbonaceous grains. Figure 
10 shows the best-fit fractional alignment /(flcff)- As previously found by Kim & Martin, silicate 
grains with b/a = \f2 must be nearly perfectly aligned for a^^ > 0.2/Ltm if the observed polarization 
is to be reproduced. 
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6. X-Ray Scattering from Partially Aligned Interstellar Grains 

6.1. Dust at a Single Distance 

We have calculated the X-ray scattering properties of graphite and silicate dust grains with the 
size distributions and fractional alignments shown in Figs. 9 and 10, for 3 X-ray energies: E = 0.5, 
1, and 2keV. The total scattering cross section Csca is given in Table 1 for the two partially- aligned 
grain models, for two viewing directions: z || B and z _L B. The two models have similar values of 
CTgca, - the total strength of X-ray scattering does not appear to discriminate between grain models 
with different degrees of grain elongation, if the models are already constrained to reproduce the 
observed optical extinction and polarization. 

Also shown in Table 1 are X-ray scattering cross sections calculated for a graphite-|-silicate 
grain model for spherical grains with the size distribution adopted by Weingartner &; Draine (2001; 
hereafter WDOl). The WDOl size distribution is based on using graphite and silicate spheres plus 
PAH molecules to reproduce the observed interstellar extinction from the infrared to the ultraviolet; 
the resulting size distribution differs in detail from what is found here for spheroids contrainted to 
also reproduce the polarization of starlight, but the total X-ray scattering cross sections cTgca are 
very similar to those obtained here for spheroidal grains. 

The ratio of X-ray scattering to visual extinction, TscaMvi is plotted in Fig. 11 for both the 
present models and the WDOl model, together with observational determinations. Aside from 
the apparent agreement between model and observation for the soft X-rays from Nova Cyg 1992 
(Draine &; Tan 2003) , there appears to be a general tendency for the observationally-determined X- 



Table 1. Total X-Ray Scattering Cross Section cxsca (lO'^"^ cm^/H)^ 
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10^^ cm^/mag 



(Bohlin et al. 1978). 
^Rv = 3.1, CpAH/H=55ppm 
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Fig. 11. — Tsca./Av VS. energy E, where Tgca is the total optical depth for scattering, and Ay is the visual extinc- 
tion. Solid curve: scattering calculated for the present size distributions for graphite spheres and silicate spheroids 
(randomly-oriented) with b/a = 1.414 (results for b/a — 2, not shown, are nearly identical). Broken curve: scattering 
calculated for graphite and silicate spheres with the WDOl size distribution. Also shown are observational results. 
For the Cyg X-1 observations by Yao et al. (2003), the filled diamonds are corrected for scattering at ^ > 120" (see 
Draine 2003b). For Nova Cyg 1992, results are shown for Tsca = 0.21 (Draine & Tan 2003) and two estimates of 
E{B - V): 0.19 mag (Mathis et al. 1995) and 0.32 mag (Vanlandingham et al. 2005). 



ray TgcaMv to be smaller than predicted by grain models based on mixtures of solid carbonaceous 
grains and solid silicate grains. However, measuring the X-ray scattering by dust is difficult at small 
angles, where it is mixed with the point spread function, and at large angies, where background 
corrections are uncertain. At this time we consider the current grain model to be viable, despite the 
fact that most of the observational data in Fig. 11 appears to fall significantly below the prediction. 
Definitive measurement of the total X-ray scattering cross section per unit Ay would be of great 
value. 

The total scattering cross sections given in Table 1 appears to decrease by ~ 6—35%, depending 
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Fig. 12. — Differential scattering cross sections per H nucleon vs. scattering angle Q for the grain mixtures of Fig. 
9 at </> = (solid line) and = 7r/2 (dotted line). Results are shown for E = 0.5, l,2keV, and for silicate grain 
axial ratio b/a = \/2 and b/a = 2. For optically-thin scattering by dust at a single distance, the halo intensity 
1(9, 4>) oc d(7/dO(0, 0), where Q = 9/(, iov (, = z/D^ where z is the distance from the source to the dust and D is the 
distance from the source to the observer. 



on energy and 6/a, when the viewing angle is changed from z ± B (maximum optical polarization) 
to z II B (zero optical polarization). This small effect would probably be difficult to confirm 
observationally. 

Fig. 12 shows the differential scattering cross section per H nucleon as a function of scattering 
angle 0, for = and = 90°, assuming the magnetic field to be in the x direction (^ = 0). For 
> 0, the scattering is stronger for = (the direction of the short axis of the aligned grains), 
but the differences between (f) = and (f) = 90° are of course reduced compared to the single-grain 
results in Fig. 2 because we are averaging over an extended size distribution, and except for the 
largest grains, there is only partial grain alignment. Nevertheless, the differences can be as large 
as ~40% at 6* > 200(keV/£;) arcsec. 
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Fig. 13. — Contours of constant cr(e, (0, 0) = 10""'"' for £ = 1 keV for grain models including partially-aligned 
silicate grains with h/a = \/2 and 2, for magnetic field B || x. 

Fig. 14 shows the anisotropy measures R^^ (0, 9) and R^^^ (0, 0) for the partially-ahgned grain 
models with b/a = \/2 and 2. Note that, for fixed 6E, the anisotropy measures R2 and R4 are 
larger at higher energies. The higher Z elements in the silicate grains (versus Z = 6 for the carbon 
grains) cause |m — 1| to decline less rapidly with energy than for carbon grains; as a result, the 
silicate grains provide an increasing fraction of the total scattering as the energy is increased. Since 
the silicate grains are aligned, but the carbon grains are not, the anisotropy measures R2 and R4 
therefore increase with increasing E. 

The top panel of Fig. 14 shows g{0,Q), the fraction of the scattered power having scattering 
angles < O, calculated for E = 0.5, 1, 2keV and the grain models for b/a = V2 and 2. If the dust 
is located at a single distance from the observer, not too close to the source, the anisotropy 
measures of the observed halo are directly related to the anisotropy measures R'f^ of the dust 
mixture: 

i^f(^„^.) = <)(|,|) (45) 

where = {D — Dd)/D, where is the distance of the dust from the observer. Thus far we have 
been considering scattering by dust at a single distance. This will apply, for example, if most of 
the dust on the line-of-sight to the source is concentrated in a single cloud and the source-cloud 
distance was large compared to the extent of the cloud along the line-of-sight. It also applies to 
scattering by Galactic dust of photons from an AGN, QSO, or GRB, in which case C ^ 1- 
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Fig. 14. — 7?<''-'(0,e) and ^^"'(0,9) vs. QE for dust mixtures with realistic size distributions and degree of 
alignment. For grains in a sheet at a distance (1 — QD from the observer, where D is the distance to the source, 
R'^\0,e) = Ri^\0,e = C©), where 9 is the observed halo angle. Results for E = 0.5, l,2keV are shown for grain 
models where the silicate grains are oblate spheroids with axial ratios b/a = \/2 and 2. 



6.2. Uniformly-Distributed Dust 

If the dust is distributed between the observer and the source, the scattering halo depends 
on the specific dust distribution. As a simple example, we consider dust with uniform density for 
0.01 < z/D < 1, with zero density for z < .OlD.^ 



^Thc dust-frcc zone near the source is introduced in order to keep /(O, 0) finite. For dust uniformly-distributed all 
the way to the source, it can be seen that I{9, 4>) oc for — > 0. This divergence is of little practical consequence, 
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Fig. 15. — Contours of constant 1(0, <p)/ I{0n,O) = 10""^^, for ^jv = 150", for uniformly-distributed mixture 

of partially-aligned grains if silicate grains have axial ratio (a) h/a — (b) hja — 2. Contours are labelled by 
I/I{9n,0). Contours plotted are for E = IkeV; for other energies, the halo intensities can be estimated from the 
approximate scaling relation d-{E, 0) w cr(lkeV, 0£/keV) 

Fig. 15 shows contours of scattered intensity for this case. In contrast to Fig. 4, the isointensity 
contours are no longer perfect elhpses for this case, as the result of the scattering contributions from 
the carbon spheres and nonaligned silicate spheroids. Nevertheless, it is seen that the isointensity 
contours are noticeably noncircular. Fig. 16 shows the normalized scattered intensity in the ^ = 
and <p = 90° directions for the two partially-aligned mixtures. 

The anisotropy measures i?2^^ (^i , 62) and R^^ {9i , 62) for uniformly-distributed partially-aligned 
dust are shown in Fig. 17, for the two values of b/a. Results are shown for a number of choices 
of 9i, the inner radius of the annulus over which 1(9, (f)) is assumed to be measured. The function 
5(^1)^2)) shown in the top panel, gives the fraction of the total scattered power that falls in the 
annulus [^1, ^2]- The observer will want to ensure that the chosen annulus [^i, ^2] will have enough 
scattered signal in it to allow reliable determination of the moment R2 and, ideally, i?4. 



because the enclosed power /o I2Tre'd9' -> as 61 -» 0. 
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Fig. 16. — I{9,4>)/I{9n,0) vs. 6 for 4> — (solid curve) and = it/2 (dotted curve) for uniformly-distributed 
grains at E = 0.5, 1, 2keV, for assumed silicate axial ratio b/a = \/2 and h/a = 2. The arbitrary normalization point 
^jv = 150arcsec/S(keV). 



7. Discussion 

7.1. Observational Considerations 

We will assume that there is no uncertainty regarding the position of the point source on 
the sky; this position is taken to be the origin of a polar coordinate system {9, ^obs) with the 
direction of (/^obs = taken to be whatever coordinate system is convenient to the observer. We 
further assume that the image I{9,<pohs) has had the instrumental point spread function (p.s.f.) 
subtracted; this is especially important if the p.s.f. is not azimuthally symmetric. Real images may 
contain extraneous background or foreground sources that need to be recognized and removed; we 
assume that standard methods are used to interpolate the scattered intensity I{6, (?!>obs) in such 
regions. In the following discussion it is assumed that I{0, (f)ohs) is the observed scattered intensity. 

An annulus [^i, ^2] is chosen that contains a strong scattered signal, with 9i chosen to avoid the 
central core that will be dominated by the instrumental p.s.f., and 62 chosen to ideally extend beyond 
the radius where the scattered halo intensity begins to significantly decline. If l{0,<pohs) is the 
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Fig. 17. — g^^\6i,9), the fraction of the scattered energy in the annulus [Si, 5], and the asymmetry parameters 

R2^\9i,8) and R^P{6i,6) vs. 9 for the annuhis [Oi.O], for grains uuiformly-distributcd between observer and source, 
for E = 0.5, l,2keV, and where the sihcate grains have been assumed to have axial ratio h/a = \/2 and 2. 



b 



observed intensity (corrected for nonuniform backgrounds), a new coordinate system (p = ^obs ~ 

A 

is defined, with the rotation angle A determined by the condition 

A = ^arctan(A/5) , (46) 

1-62 i-2tt 

A = / dOsine #obs^(^,</'obs)sin(20obs) , (47) 
Jei Jo 

B = / dOsmO #obsi"(^, <^obs) cos(20obs) • (48) 
JOi Jo 

(49) 

In this coordinate system {0,(j)), the anisotropy statistics R^p {61,62) can be calculated from eq. 
(36). The resulting (61,62) and R^ {61,62) can be compared to the theoretical predictions in 
Figure 14 (for dust at a single distance) or Figure 17 (for uniformly-distributed dust). 

We have calculated the anisotropy statistics R^ for plane-parallel distributions of dust: dust 
in a sheet at a single distance and dust uniformly distributed between observer and source. In real 
situations, gradients in the dust density distribution perpendicular to the line-of-sight (on linear 
scales of order ^haio ~ 1' x Ikpc ^ 0.3 pc) can contribute to the R^ ■ How can halo anisotropies 
due to aligned dust grains be separated from anisotropies due to dust density gradients? 
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1. The halo anisotropy due to dust grain ahgnment wih have = R3 =0. Dust density 
gradients, on the other hand, wih generally contribute to all of the R[^\ and in fact would 
be expected to make the largest contribution to r[^^ . Therefore the magnitude of R^^ and 
i?3^^ can be used to estimate the contribution of dust density gradients to .^2^^ 

2. The halo anisotropy due to aligned dust grains will be aligned with the direction of starlight 
polarization: polarized starlight has E parallel to the short axis of the aligned grains. There- 
fore the scattered halo will have its major axis parallel to the direction of starlight polarization: 
the angle A should coincide with the direction of starlight polarization. 

3. If the anisotropic signal i?2^^ is from aligned grains, the angle A will be insensitive to the 
choice of annulus [^1, ^2]- If A is found to vary significantly from one annulus to another, the 
observed anisotropy will have an appreciable component from some source other than aligned 
dust grains. - nonuniform dust, or other X-ray sources in the field. 



Ideally, the X-ray scattering halo would be observed for a source where the optical starlight 
polarization is known cither for the source itself or a nearby stellar companion. In some cases, the 
X-ray source may not be bright enough at optical wavelengths to permit polarization measurements. 
However, the polarization of starlight is known to display large-scale coherence over the sky (see, 
e.g., Mathewson &; Ford 1970); therefore if starlight polarization has been observed for a star that 
is nearby on the sky, with a similar amount of reddening, it is reasonable to presume that the dust 
in front of the candidate X-ray source would produce a similar degree of optical polarization. 

Directions with large observed starlight polarization tend to lie close to the galactic plane 
(so that significant amounts of reddening will be present) and in directions perpendicular to the 
general direction of the local magnetic field. An example of a region with large observed starlight 
polarization per unit E{B — V) is I = 130 it lOdeg. Figure 18 shows piV) / E{B — V) for stars in a 
section of the Galactic plane near I « 130 deg. The stars shown are limited to 0.30 < E{B — V)< 
0.60 mag (i.e., 0.9 < Ay 1-9 mag). The large-scale coherence of the polarization direction is 
apparent, and it is also evident that the ratio of polarization to reddening is relatively uniform. 
This indicates that the magnetic field is relatively well-ordered, and that the efficiency of grain 
alignment is also relatively uniform. An X-ray source in this general region would presumably be 
observed through the same aligned dust as the stars for which the polarization has already been 
measured. 



7.2. Predictions 

We have considered two models for interstellar grains, both with partially aligned silicate 
oblate spheroids, but for two different axial ratios. The aligned grains that produce polarization 
of starlight also produce anisotropic X-ray scattering. While we have explicitly calculated the 
anisotropic X-ray scattering only for the optimal case where the magnetic field B is unidirectional 




Fig. 18. — Line segments show ratio direction of starlight polarization, with length proportional to p{V)/E{B — V), 

where p(V) is the polarization, and E{B — V) is the reddening. Data from Hciles (2000). In this region of the sky, 
the dust grain ahgnment is evidently quite uniform, and the polarization per unit E{B — V) is ~ 2/3 of the empirical 
maximum. 



and perpendicular to the line-of-sight, wc expect that for other magnetic field geometries, the 
anisotropy i?2^^ will decrease in proportion to Pma.x/E{B — V), because for optically-thin conditions 
both i?2^^ and Pmax/E{B — V) respond similarly to changes in the magnetic field direction.'^ We 
can therefore predict the ratio of the fractional anisotropy of the X-ray halo to the ratio of the 
optical polarization to total reddening. 

We suppose that it is possible to measure the X-ray halo for 9i < 6 < 82, with 61 < 200", and 
02 ^ 300". If the dust is distributed approximately uniformly between observer and source, from 



'^E.g., both are unchanged if B — » — B; both vanish if B || z; and both vanish if 50% of the dust has symmetry 
axis ai || x, and 50% has ai || y. 
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Figure 17 we see that for X-ray energy E = 1 keV, 

^ , (50) 



4'^(^i,^2)>0.09x 



OmE{B -V)/mag 



provided the annulus has 9i < 300", and 62 > 300". The two cases b/a = \/2 and 2 give similar 
values for the quadrupole anisotropy i?2^^ - increased grain oblateness is offset by reduced degree 
of partial alignment. As discussed above, the anisotropy r'^'^ is an increasing function of energy 
over the range 0.5-2 kcV because the relative contribution of the nonspherical silicate grains to the 
scattering increases with energy. 

It is not critical that the dust be distributed uniformly for the result (50) to be applicable, 
only that a negligible amount of the starlight polarization be produced by dust that is very close to 
the source, as the anisotropic X-ray halo produced by this dust may be largely lost at small angles 
unless it is possible to use small values of Oi. 

Aligned grains are also predicted to produce an octupole anisotropy in the X-ray halo, mea- 
sured by Unfortunately, this anisotropy can be washed out by moderate rotation of the 
magnetic field direction along the line-of-sight. However, if the starlight polarization has Pmax ~ 
0.09E{B — F)/mag, then we can suppose that the magnetic field in the dusty regions is approxi- 
mately perpendicular to the line-of-sight and unidirectional (otherwise the net starlight polarization 
would be reduced). In this case, we can use the statistic R^^^ as a test of grain models. The model 
with minimal axial ratio b/a ^ 1.4 for the spheroids predicts R^/^ « 0.03 at 1 keV (see Figure 
17). If larger values of R^/^ are observed, it will suggest more extreme grain shapes (for example, 
R^^ 0.10 would be consistent with b/a « 2 for the silicate spheroids, but inconsistent with 
b/a ^ ^2). Therefore, measurement of R^^ has the potential to discriminate between grain models 
with different grain shape. 



8. Summary 

The principal results of this paper are as follows: 

1. We show that anomalous diffraction theory can be used for accurate calculation of X-ray 
scattering from dust grains with arbitrary shape. Fast Fourier transforms may be employed 
for efficient calculation of the scattering halo. 

2. Differential scattering cross sections are calculated for oblate spheroids, and are shown to 
have substantial anisotropics if the spheroid axial ratio differs appreciably from unity. 

3. Models of spherical carbonaceous grains and partially-aligned silicate spheroids are found 
that reproduce the observed interstellar extinction and polarization of starlight as a function 
of wavelength. Confirming a result found previously by Kim & Martin (1995), we find that 
for silicate axial ratio b/a \/2 the model can reproduce the largest observed starlight 
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polarization, provided that the sihcate spheroids with a^s > 0.1/xm are almost perfectly 
aligned (with their short axis parallel to the magnetic field). For larger axial ratios, the grain 
alignment need not be so complete. 

4. X-ray scattering halos are calculated for aligned interstellar dust grains. For realistic size 
distributions and fractional alignments, the scattered halo shows substantial anisotropy. 

5. We propose statistics R^P {61,02) to measure the halo anistropy in an annulus 9i < 6 < 62- 
We predict the values of R^^ /\pinax/E{B — V)]. We find that R^ is large enough to be 
measured on sightlines to X-ray point sources where the grains are aligned so as to produce 
starlight polarization with Pmax/E{B — V)> 0.05/mag. 

6. We show that the octupole anisotropy R^ is sensitive to the assumed grain shape. On 
favorable sightlines, R^ can be used to constrain the geometry of interstellar grains. 

This research was supported in part by NSF grants AST-0216105 and AST-0406883. We are 
grateful to an anonymous referee for helpful comments. B.T.D. is grateful to R. Bandiera for 
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A. Least-squares fitting procedure 



Let nc(a) be the number density of grains with composition c and radii < a. It is convenient 
to define u = Ina and yc such that dric/dlna = dric/du = exp(yc). Then duc/da is positive definite 
for —GO < Uc < oo. 

The size distribution is sampled at N^ad sizes, and extinction and polarization are calculated 
at Nx wavelengths. 



The terms include the following contributions 
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The radii aj arc assumed to be uniformly distributed in u = In a: Uj^i — uj = Au. The deriva- 
tives are evaluated using the usual differencing: {df /du)j^i/2 = {fj+i — fj)/^u and {d^f/du?)j = 
(/,•+! + /,_i-2/,)/(A«)2. 

Here KiiLO = 2.29 x lO'^^ cm^/H and l^,ar,o = 1-57 x lO'^^ cm^/H are target values based on the 
estimated elemental abundances of species such as Mg, Si, Fe, and C, the fractions of these elements 
believed to be in solid form, and assumed densities for the silicate and carbonaceous material. For 
the target volume for silicate grains we assume a composition Mgi.iFeo.9Si04, with Mg/H=3.3 x 
10"^, and density p = S.Tgcm"^. This would consume 87%, 93%, 94%, and 29% of the current 
estimated solar composition values for Mg, Fe, Si, and O, respectively [(Mg/H)0 = 3.8 x 10^"'' 
(Grevesse & Sauval 1998), (Fe/H)0 = 2.9 x IQ-^ (Asplund et al. 2000), (Si/H)0 = 3.2 x lO^^ 
(Asplund 2000), and (O/H)© = 4.6 x 10"*^ (Asplund et al. 2004)]. For the target volume of 



-35- 



carbonaceous grains we assume 70% of the current estimated solar abundance (C/H)0 = 2.46 x 10 
(Asplund et al. 2005) to be in grains with p = 2.2gcm~^. 

The coefHcients ai and a2 weight the penalties for deviating from the target abundances; 
as and a4 weight the penalties for fractional errors in reproducing the observed extinction and 
polarization; weights the penalty if the alignment fraction f{a) is not a monotonically- increasing 
function of size a; uq and cty penalize non-smoothness in the size distributions; and ag penalizes 
non-smoothness in the alignment function. 



